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We study phase transitions in homogeneous spin-1 Bose gases in the presence of long-range magnetic 
dipole-dipole interactions (DDI). We concentrate on three-dimensional geometries and employ momentum 
shell renormalization group to study the possible instabilities caused by the dipole-dipole interaction. At the 
zero-temperature limit where quantum fluctuations prevail, we find the phase diagram to be unaffected by the 
dipole-dipole interaction. When the thermal fluctuations dominate, polar and ferromagnetic condensates with 
DDI become unstable and we discuss this crossover in detail. On the other hand, the spin-singlet condensate 
remains stable in the presence of DDI. 
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I. INTRODUCTION 

The past few years have revealed that ultracold atomic gases 
can answer important questions beyond the immediate scope 
of atomic physics Dffil. In particular, experimental methods 
have matured to the level where measurements of critical ex- 
ponents are possible in some cases |2|. This provides an in- 
teresting opportunity to study the physics of phase transitions 
and critical phenomena utilizing cold atomic gases, as well as 
to realize exotic phases that are absent in more conventional 
solid state systems 0. In this work, we consider Bose gases 
with a spin degree of freedom (SI |7) . They provide an in- 
triguing example where magnetic ordering can compete with 
superfluidity and condensation. This interplay can give rise to 
a myriad of topological defects llBT flBI which play an impor- 
tant role, e.g., in the superfluid transition in low-dimensional 
systems lfT71[T8l . 

Initially, the magnetic properties of spinor Bose gases were 
assumed to depend only on the local interactions determined 
by the scattering lengths in the different total hyperfine spin 
channels fl6j[7jD9-- 22 ], but recent experiments suggest that the 
long-range magnetic dipole-dipole interaction (DDI) may be 
an essential ingredient in determining the properties of spinor 
Bose gases l23Tl27ll . In this work, we consider the effect of 
dipole-dipole interaction in spin-1 Bose gases using momen- 
tum shell renormalization group (RG) [28-33 1. We note that 
also the functional renormalization group [34] has been suc- 
cessfully applied in the context of cold atoms 0514371 . The 
momentum shell RG analysis allows us to determine the ef- 
fect of DDI on the phase diagram of spin-1 Bose gases which 
has recently attracted some interest [38-40]. Moreover, the 
recent advances in the creation of Feshbach resonances using 
either optical means [41 ] or microwaves [42 1, suggest that ex- 
ploration of the phase diagram could become experimentally 
realistic in the near future. 

Dipole-dipole interaction couples the spin directly to spa- 
tial degrees of freedoms, giving local spins tendency to align 
head to tail and antialign side by side ||43~1 l44l . On the other 
hand, the experiments described in Refs. [26, 27] are of mixed 
dimensionality in the sense that spin dynamics was effec- 



tively two-dimensional while otherwise the system was spa- 
tially three-dimensional (3D). Furthermore, the original DDI 
was strongly modified by a rapid Larmor precession induced 
by an external magnetic field. In the present work, we fo- 
cus on the properties of pristine DDI and consider a homo- 
geneous three-dimensional system in the absence of external 
magnetic fields. For three-dimensional systems, DDI is a true 
long-range interaction [44, 45] and we avoid additional com- 
plications that may arise due to the absence of true long-range 
order in low-dimensional systems. 

Although dipole-dipole interactions are present in all ferro- 
magnetic materials, they are usually weak and often neglected 
or treated phenomenologically ll46l H71 . However, for fer- 
romagnets which order only at very low temperatures, DDI 
might be crucial for the correct low-energy behavior 1481 . In 
this work, we analyze this scenario in the context of spin-1 
Bose gases. We find that DDI introduces additional instabili- 
ties to the expected finite-temperature phase diagram 11381 . In 
particular, DDI renders both polar and ferromagnetic conden- 
sates unstable and the RG analysis alludes to the existence of 
a fluctuation-induced first-order transition. 

In the zero-temperature limit, we show that DDI renormal- 
izes to zero and the usual mean-field theory |6] [7] |44) is a 
valid description of the system. Dipole-dipole interactions 
also generate a new single-particle term which has not been 
taken into account in the previous studies. This new interac- 
tion is relevant in the RG sense and it is allowed by the sym- 
metries of the system. However, in the zero-temperature limit 
it renormalizes to zero along with the DDI. 

II. THE MODEL 

We consider a uniform spin-1 Bose gas neglecting the ef- 
fects of an external potential that confines the atoms. In the 
presence of DDI, the system has a global U(l) symmetry as- 
sociated with the conserved atom number and a global SO (3) 
symmetry corresponding to a simultaneous rotation of spin 
and coordinate spaces |49l |50l . The latter symmetry indi- 
cates that only the sum of spin and orbital angular momen- 
tum is conserved. The effective action in the Zeeman basis 
{\F = l,m F = +1,0, -1)} can be written as S = S Q + S int , 
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where r is the imaginary time and j3 = l/k B T. We always 
assume implicit summation over repeated indices. The local 
spin is given by 5* = F ab ip b , where T % are spin-1 ma- 
trices in the Zeeman basis and ip a , a = —1, 0, 1 are bosonic 
fields. Parameter V is initially set to unity and it acquires non- 
trivial renormalization under the RG transformation. In this 
work we consider two distinct limits: the T = case where 
r renormalizes only due to the anomalous dimension of the 
fields tp a and the high temperature limit where T renormal- 
izes to zero and we obtain a classical theory. 

The bare values of the coupling constants c and c 2 are re- 
lated to scattering lengths a Q and a 2 in the total hyperfine spin 
channels F = and F = 2 by c = 4nh 2 (a + 2a 2 )/3m 
and c 2 = 4irh 2 (a 2 — a )/3m. The coupling constant corre- 
sponding to DDI is given by c dd = fi fi 2 g 2 /4TT, where fj, is 
the vacuum permeability, /i B Bohr magneton, and g F Lande 
g-factor. The kernel for dipole-dipole interactions in the mo- 
mentum space takes the form [43 1 
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In real space, ^(r) decays as l/|r| 3 . To date, the only ex- 
perimentally studied dipolar spin-1 Bose gas has been 87 Rb 
for which the different coupling constants satisfy c 2 /c = 
—0.005 and c dd /|c 2 | = 0.1 G6l . Another candidate for dipo- 
lar spin-1 Bose gas is 23 Na for which the scattering lengths of 
Ref. ||5T| give c 2 /c = 0.03 and c dd /c 2 = 0.006. The small 
value of c dd /c 2 explains why the effects arising from DDI 
have been expected to be vanishingly small for 23 Na. The dif- 
ferent interaction vertices appearing in the RG calculations in 
Sections III VI are illustrated in Fig. [T] 

To streamline the RG calculations, we switch to the Carte- 
sian basis. In this basis, the field operator ^ = (ip x , ipy,ip z ) 
transforms as a vector under spin rotations. Moreover, the 
spin- 1 matrices ( T x ,F y ,J- z ) take a particularly simple form 
(J r Q ) ) g 7 = — ie a pj, where e Q) g 7 is the Levi-Civita tensor. The 
effective action can be written as 
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where indices {x, y, z} are referred to by the Greek indices 
a, (3, ... and the Latin indices a,b, ... correspond to the origi- 
nal Zeeman basis. We use a shorthand notation k = (u; n , fe) 
and / d 4 fc = J dk. Bosonic Matsubara frequencies are 
given by u) n = 2irn/h(3 and e k = h 2 k 2 /2m. The coupling 
constants g d and g s are related to the coupling constants in the 
Zeeman basis by g d = c + c 2 and g s = — c 2 . 



III. RENORMALIZATION GROUP CALCULATION 

We set up the RG calculation in a fixed dimension D fol- 
lowing Refs. Il30l433l . To make a connection to Refs. l38l[39l , 
we first neglect the dipole-dipole interactions and study a gen- 
eral £>-dimensional situation. We show that our RG equations 
at the zero-temperature limit coincide with Ref. 1 39 ] and es- 
sentially reproduce the phase diagram proposed in Ref. Il38l . 
We also point out that in a contrast to Ref. Il39l . in which 
the stability of low-dimensional multi-component Bose gases 
was considered at the zero-temperature limit, our main fo- 
cus is a three dimensional spinor Bose gas at finite temper- 
atures. We note that isotropic long-range interactions in spin- 
less Bose gases have been analyzed in the zero-temperature 
limit in Ref. l32l and long-range interactions of the form 
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Figure 1. (a) Generic interaction vertex, (b) local density-density 
interaction, (c) local spin-spin interaction, (d) and dipole-dipole in- 
teraction. The generic interaction vertex (a) can denote any of the 
vertices (b)-(d). Conservation of the 4-momentum k = (oj„, k) at 
each vertex is implied. 
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V(r) oc 1/H S were found irrelevant for s > 2. Our find- 
ings in the presence of DDI are similar to those of Ref. 11321 . 
namely, DDI becomes irrelevant at zero temperature. 

To study the effects of DDI, we employ the momentum 
shell RG 1 28 . 29] EH m which we split the fields appearing 
in Eqs. Q and §5§ such that ip a = if; a < + ip a >, where the 
ip a < contains momentum components with |fe| < A/s and 
tp a > corresponds to momenta A/s < |fc| < A. The ultravio- 
let (UV) cutoff is denoted by A, and in general, nonuniversal 
quantities such as condensate fraction or critical temperature 
depend explicitly on A. Several tricks such as halting the RG 
flow when an appropriate scale is reached or relating A to the 
s-wave scattering length can be used to obtain information on 
quantities depending on A 130 -331 1391 . 

The RG calculation proceeds by integrating out the fast 
modes residing at the momentum shell A/s < \k\ < A, which 
results in a renormalized action for the slow modes ip a < . 
At the second step of RG transformation, the UV cutoff is 
brought back to the original value A by rescaling the fields 
and momenta, giving rise to RG equations for the chemical 
potential and coupling constants. Only the one-particle irre- 
ducible connected diagrams contribute to the RG equations. 
In this work, we compute the RG equations to the one-loop 
order. The relevant diagrams appearing in the renormaliza- 
tion of chemical potential and coupling constants are shown 
Figs.|2]and|3] 
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Figure 2. One-loop diagrams contributing to renormalization of the 
chemical potential. The interaction line can be any of those denoted 
in Fig. l(b)-(d). The external legs corresponding to fields ifr a < and 
ipa,< are denoted by gray lines for clarity. The dipole-dipole interac- 
tion can give rise to external interaction lines which are also denoted 
by gray dashed lines [the tadpole diagram in (a)]. 



where we have taken s — e e . For simplicity, we first ne- 
glect the anomalous dimension of the fields, which allows us 
to keep the kinetic energy term in Eq. Q fixed during the RG 
transformation. In Section [VI] we take into account also the 
renormalization of the kinetic energy term and find that e k is 
only weakly renormalized. For vanishing anomalous dimen- 
sion we obtain an identity 1 30 1 



2C + z = 2 - D, 



(8) 



and the relevance of all other terms is compared to the kinetic 
energy. The requirement that the rescaled action is equivalent 
to the original one yields the scaling relations 
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where g = g d , g s , c dd , and we have used Eq. dHJ. In the pres- 
ence of DDI, we take D = 3. 
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Figure 3. One-loop diagrams contributing to renormalization of cou- 
pling constants g d , g B , and c dd . The two interaction lines can be any 
combination of the interactions in Fig. l(b)-(d). The external legs 
and interaction lines are again denoted by gray color. 



The diagrams in Figs. [2] and [3] correspond to an expansion 
with respect to coupling constants g d , g s , and c dd in the first 
non-trivial order. The internal lines are evaluated using non- 
interacting one-particle Green's function 
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After integrating out the fast modes and neglecting the ir- 
relevant terms generated by the momentum shell integration, 
the slow fields, momentum, and imaginary time are rescaled 
as [30| 
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In general, RG calculations provide information on univer- 
sal quantities such as different phases and transitions between 
them. On the other hand, renormalization group analysis can 
be used to determine the relevance of a particular interaction 
for a given phase and to study the stability of different phases 
when additional interactions are included. We take the latter 



(6) point of view in Sections IV VI where we analyze the effects 



of DDI on the phase diagram of spin-1 Bose gases. 



IV. RG FLOW IN THE ABSENCE OF DIPOLE-DIPOLE 
INTERACTIONS 

Let us first consider a general dimension D and neglect 
DDI. The different diagrams contributing to the RG equations 
are shown in Fig. |4] Since all interactions are local, the dif- 
ferent diagrams in Fig. [4] reduce to evaluation of the bubbles 
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shown in Fig. [5] Integration on the momentum shell is re- 
stricted to the interval A/s < |fc| < A and the Matsubara 
sums can be calculated using the standard methods |52"1 . At 
the limit of an infinitesimal shell of thickness AM, we obtain 

dfi = -2{g s + 2g d )K D T n BE [f3T(e A - fi)]A D d£, (10) 
d.9 s = -[(3 5s 2 + 2g d g s ) Xl + 4 gd9sX2 ]K D A D de, (11) 
d<? d = ~[9dXi + 4(5d + 9 d 9 s + 9s)x 2 }K D A D de, (12) 

where K D = [2 D - l n D / 2 T(D/2)\- 1 and T(x) is the Gamma 
function (not to be confused with the energy parameter T). 
Furthermore, functions Xiifti T) and X2(/3> T) are given by 
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Figure 5. Basic bubbles for diagrams in the absence of dipole-dipole 
interactions. The bubbles depend in principle on external momenta. 
Expanding about the zero momentum one observes that terms cor- 
responding to finite external momentum are at most marginal in the 
RG sense. Since we have neglected marginal interactions from the 
outset, we take the external momenta to be zero. 



xi(/3,r) = r 
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equations 
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We have denoted the Bose distribution function by n SE {x) = 
l/(e x — 1) and e A — h 2 A 2 /2m. The contributions propor- 
tional to Xi correspond to diagrams containing the bubble in 
Fig.|5ja) and contributions containing Xi arise from diagrams 
with the bubble in Fig.|5jb). 







Figure 4. One-loop contributions to RG equations in the absence of 
dipole-dipole interactions. 



Taking into account the scalings dictated by Eqs. (|7]i 
and using the scaling relation in Eq. (jH), and transform- 
ing back to the Zeeman basis we obtain the following flow 
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-zf3, (15) 

-(2( + D)T, (16) 

2 M - 2(c 2 + 2c )^A Z3 rn BE [/3r( £A - fj,)], (17) 

= 2(C + l)c - [(eg + 2cl) Xl + 4c 2 oX 2}K D A D , (18) 

= 2(C + l)c 2 - [(2c c 2 - c 2 2 ) Xl 

+ A(c Q c 2 +c 2 2 )x 2 }K D A D . (19) 



Although these equation are valid for any temperature, we 
consider two limits: the quantum regime which takes place at 
the zero-temperature limit and is dominated by the quantum 
fluctuations, and the thermal regime where thermal fluctua- 
tions prevail over the quantum fluctuations 1301 . 



A. Quantum regime 

Let us first consider the quantum regime at which we re- 
quire dT/d£ — 0. Furthermore, we set V = 1. From Eqs. ^ 
and ( fl6] l we obtain £ = —D/2 and z = 2. This gives the 
usual instability of the T* = fixed point, since any nonzero 
temperature tends to increase in the RG flow. In the limit 
k B T e A — /i, we have 



xi(/3,r) 



2(£ A - M )' 

xa(&r)->o, 

where the latter equation gives the well-known result stating 
that in the zero-temperature limit, only the ladder diagrams 
contribute to the renormalization l53l . The RG equation for 
the chemical potential becomes d^i/dl = 2/z, and the only 
fixed point is n* — 0. Furthermore, the remaining RG equa- 
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dc o = 
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dc„ 
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K D A D 
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(2-D)c 2 + (c 2 2 -2c c 2 ) K f A ^. (21) 
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The above equations are precisely those of Ref. ||39l , and for 
the future reference, we consider here the D — 3 case. Cases 
D = 1 and D = 2 have been analysed in Ref. 
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Table I. Different fixed points corresponding to the RG equa- 
tions {20\ and {21) . The dimensionless values are defined as c* = 
2£aC* / K D A D , i = 0, 2. In three dimensions, the Gaussian fixed 
point I is stable and the SU(3) symmetric fixed point IV is unstable. 
Fixed points II and III have both relevant and irrelevant scaling fields. 

The fixed points corresponding to the RG flow defined by 
Eqs. ( |20] i and pTj ) can be determined exactly. The four dif- 
ferent fixed points are given in Table [I] and the RG flow for 
D = 3 is shown in Fig. [6] The fixed points Cq and c 2 corre- 
spond to the dimensionless values c* — 2e A c* /K D A D , for 
i = 0, 2. In the special case D = 2, the dimensionless quan- 
tities are independent of A [39]. The RG flows in Fig.[6]show 
that similarly to the D = 1 and D = 2 cases, there are two 
runaway flows indicating the formation of bound spin singlet 
pairs (positive c 2 ) and ferromagnetic instability (negative c 2 ) 
where the condensate becomes locally fully spin-polarized in 
the sense that fluctuations in the magnitude of the local spin 
become suppressed. 

The runaway flow associated with the formation of pair 
condensate renders the coupling g = c Q — 2c 2 correspond- 
ing to the spin singlet channel negative, while the coupling 
g 2 = c + c 2 in the F = 2 channel remains positive. Hence 
this instability corresponds to formation of spin singlet pairs. 
The second runaway flow where c 2 becomes ever more nega- 
tive renders both g and g 2 negative. We refer to these two 
instabilities as antiferromagnetic (AFM) and ferromagnetic 
(FM) runaway flow, respectively. 

At mean-field level, stability of a finite cloud against col- 
lapse requires the bare coupling constants to satisfy c > 
and — c 2 < c ll38l . On the other hand, the flow diagram in 
Fig. ^indicates that many-body corrections yield a larger win- 
dow of coupling constants c and c 2 which renormalize to the 
Gaussian fixed point. Although both FM and AFM runaway 
flows suggest that the system becomes unstable, stability can 
be restored by including the higher order terms generated by 
the RG transformation. Such terms in general are marginal or 
irrelevant in the RG sense, but can nevertheless become im- 
portant in the regime where RG flow does not converge to any 
fixed point 15411551 

Between the regions corresponding to AFM and FM run- 
away flows, the gas forms the usual spinor condensate. Since 
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Figure 6. (Color online) Renormalization group flow for a three 
dimensional spin-1 Bose gas at zero temperature in the absence of 
dipole-dipole interactions. The dimensionless coupling constants are 
defined as c i — 2e A c i /K- i k i , i — 0, 2. The Gaussian fixed point (I) 
is stable, and RG flows starting from c > and small |c 2 1 end up to 
the Gaussian fixed point. Fixed points II and III are in the physically 
inaccessible region but they still have an effect on the RG flow since 
they give rise to two runaway flows for large | ti \ . 



interactions tend to renormalize to the Gaussian fixed point, 
the Bogoliubov mean-field theory of spinor condensates is a 
valid description of the system in this regime (cf. Ref. Poll ). 
However, the RG approach used here does not provide in- 
formation about nonuniversal properties such as the possible 
fragmentation of the condensate for antiferromagnetic inter- 
actions 1561. 



B. Thermal regime 

In the thermal regime, we require that /3 in Eq. ( fT3] > does 
not flow, i.e., the temperature is kept constant in the RG flow. 
This implies z = 0, and Eq. ([H) gives £ = (2 — D)/2. From 
Eq. ( fTo} we observe that any finite initial T flows to zero. 
Quantum fluctuations are thus negligible in this limit and we 
take T ^ in Eqs. <[17HI9), We obtain 
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and at the critical plane corresponding to fi = f57] [58J we 
have 



dl 
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At finite temperatures we define dimensionless coupling con- 



stants by Cj = fielcJKpA 1 



0, 2. Note that the defini- 
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tion is slightly different from the zero-temperature case. Fixed 
points corresponding to the dimensionless coupling constants 
are shown in Table [II] where have defined e = 4 — D. 
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Table II. Fixed points corresponding to the RG equations {25) 
and (26) . The dimensionless values are defined as c* = 
f3e 2 A c* / K D K D , i — 0,2. The Gaussian fixed point I is unstable 
and the SU (3)-symmetric fixed point IV is stable for e > {D < 4). 
Fixed points II and III have both relevant and irrelevant scaling fields. 

To analyze the RG equations in the absence of DDI, we 
concentrate on the case D = 3 for which the RG flows are 
depicted in Fig. [JJ The instabilities indicated by the run- 
away flows have the same structure and physical interpretation 
as in the zero-temperature case. An interesting difference to 
the zero-temperature limit is that the ferromagnetic instability 
corresponding to the runaway flows for large negative c 2 takes 
place before the mean-field criterion c > and — c 2 < c 
is violated. Hence, the thermal fluctuations tend to decrease 
the stability of spinor condensates on the ferromagnetic side 
(c 2 < 0). 

We note that only the ratio c 2 /c is universal (i.e., indepen- 
dent of the cutoff A), and therefore quantitative comparison of 
the singlet condensate formation in the quantum and thermal 
regimes depends on A, see Figs.[6]and[7j The values c 2 /c dis- 
cussed in Sec.|n]place the bare coupling constants c (0) and 
c 2 (0) for 23 Na and 87 Rb into the regime where c (£) and c 2 (£) 
either renormalize to zero (quantum limit) or to the SU (3) 
symmetric fixed point where c 2 = (thermal limit). Since 
ferromagnetic and antiferromagnetic (polar) condensates cor- 
respond to different symmetries, they should be separated by a 
phase transition. At mean-field level, we expect the transition 
to be first order (see also Ref. [38]). The RG calculation sup- 
ports this conclusion in the sense that we do not find a critical 
point separating the two phases, see Fig. [JJ 

At first it may seem surprising that the fixed point IV (or 
the Gaussian fixed point I at T = 0) governs the properties of 
both antiferromagnetic and ferromagnetic condensates. How- 
ever, according to the exact theorem of Ref. [59], the correct 
low-energy theory of spin-1 bosons should indeed correspond 
to c 2 = 0. The authors of Ref. Il59l propose that a nonzero 
c 2 in the low-energy theory could arise either from dipole- 
dipole interactions or from the electron transfer between the 
atoms. In Sees. [V] and [VI] we show that DDI is indeed suffi- 
cient to give rise to nonzero c 2 in the low-energy description 
of spinor Bose gases. Since several important properties of the 
spinor Bose gases hinge upon the presupposition that the spin- 
dependent coupling is nonzero l8T [T5l[T7l[T8l . the analysis in 
Sees . [V] and | VI|provides justification for this key assumption. 
The two runaway flows do not contradict the aforementioned 
theorem since they correspond to the formation of a spinless 
condensate consisting of either single atoms polarized to the 
same hyperfine spin state or spin singlet pairs. 
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Figure 7. (Color online) Renormalization group flow for a three di- 
mensional spin-1 Bose gas at fixed temperature. The RG flows corre- 
spond to the case where dipole-dipole interactions are either absent 
or take the fixed-point value c*[ d = 0. The dimensionless coupling 
constants are defined as c i — fie] v c i j A" 3 A 3 , i = 0,2. The ferro- 
magnetic condensate corresponds to the region with c 2 < above 
the lower dashed line. The region of antiferromagnetic condensate is 
delimited by c 2 > and the upper dashed line. 



V. RG ANALYSIS OF DIPOLAR BOSE GAS 

We calculate contributions arising from DDI using the 
generic diagrams depicted in Figs. [2] and [3] The contribution 
from the tadpole diagram in Fig. |2ja) vanishes identically for 
DDI, and the rainbow diagram in Fig.[2|b) does not contribute 
to the renormalization of chemical potential. However, the 
rainbow diagram does in principle contribute to the anomalous 
dimension rj which indicates the importance of renormaliza- 
tion of the kinetic energy term in Eq. At this point, we 
neglect the renormalization e k — > Z v e k altogether and set 
Z„ = 1. We will justify this assumption after we have derived 
the full RG equations in Sec. IV. 

The rainbow diagram also generates a relevant term which 
in the Cartesian basis takes the simple form 

S'q = -h J d 4 k ip* a {k) k a k p i/j (k), a,fi e {x, y, z}. 

(27) 

We have introduced here a new coupling constant h Q that is 
determined by the RG equations. At the lowest order, h is 
proportional to c dd . We note that the combined contribution 
of the kinetic energy renormalization and S' a can be written as 
a squared spin-orbit interaction 

S so = / d 4 krjk)(k-J^)l^ (k). (28) 

The new single-particle term S'q is allowed since it has the 
same symmetries as the original action given by Eqs. Q 
and (151. Whether this term should be included to the original 
action from the beginning depends on the behaviour of c dd (£) 
in the RG flow and, in particular, on the values of c dd (£) at the 
fixed points. 
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an addition that c 





Figure 8. Additional one-loop diagrams giving a nonzero contribu- 
tion to the RG equations in the presence of dipole-dipole interac- 
tions. 



Most of the diagrams in Fig. [3] do not contribute to the 
renormalization of the coupling constants. All the relevant 
terms are shown in Fig. [8] and they can be evaluated using 
essentially the same methods as in Sec. [TV] They yield 



d.9 s = (AidjXi + A 2 c 2 d /3x 2 )X 3 A 3 d^ 
d.9 d = ~(A 3 4lXi + \ 4 c 2 dd (3 X2 )K 3 A 3 de 



(29) 
(30) 



dc dd = -(A 5 c dd + 2c dd g d - 4c dd g s )K 3 x 2 A M, (31) 



= 272tt 2 /45, A 3 = 256tt 2 /45, 
8ir/3. The complete set of RG 



where A x = 32tt 2 /45, A 2 
A 4 = 496vr 2 /45, and A 5 = 
equations for a dipolar Bose gas consists of contributions from 
Eqs. (p9]l-(|3T| and the RG equations (p~5|>— (fT9]> . We consider 
again the same limits as in Sec. IV namely, the quantum limit 
T -> and the thermal limit T -> 0. 



A. Quantum regime 



Using the results of Sec. |IV A| we obtain the RG equations 
for nonzero dipole-dipole interaction at the fixed point /i* = 



dco 
d£ 
dc 9 



K A 3 

-c -( C 2 + 2c 2 + 7l c 2 d )- 3 



2s 



-c 2 + (c 2 - 2c c 2 - 7 2 c^ d ) 



A 

2e A : 



dc. 



(id 



b ddJ 



(32) 
(33) 
(34) 



where 7l = 2247r 2 /45 and 7 2 = 327r 2 /45. We observe im- 
mediately that the dipole-dipole coupling constant c dd renor- 
malizes exponentially to zero in the zero-temperature limit, 
and the new single-particle term Eq. ( |27] i remains unimpor- 
tant. The fixed points of the RG flow are shown in Table|I]with 



da 



at each fixed point. Furthermore, 
since c dd appears quadratically in Eqs. ( |32| i and ( [33] l, the sta- 
bility of the fixed points is not affected by DDI and c dd gives 
rise to an irrelevant scaling field at these fixed points. Other- 



wise the properties of RG flows are the same as in Sec. IV A 



B. Thermal regime 

In the thermal regime, we require again that the temperature 
does not flow under RG, which renders the parameter T to 
renormalize exponentially to zero. Using Eqs. (|22]i-(|24|i, we 
obtain, for fi* = 



= c - (5c 2 , + 2c 2 + a c 2 d ) 



/3e 2 



dc 2 
dl 
dc. 



(3c 2 . + 6c c 2 + a 2 c 2 dd ) 



A 

^A 3 



dd 



d(: 



c dd ~ ( a dd c dd + 2c c dd 



6c 2 c dd ) 



^ 3 A 3 

PA ' 



304tt 2 /45, and a 



(35) 



(36) 



(37) 



8tt/3. 



where a = 4487r 2 /45, a 
The RG flow gives rise to the four fixed points discussed in 
Sec. IV B corresponding to c dd = 0. There are also two addi- 



tional fixed points given by 



c* a = 0.087, c* 2 = 0.117, c* dd = 0.015, 
c* Q = 0.088, c 2 = 0.148, c* dd = -0.008. 



(38) 
(39) 



The dipole-dipole interaction introduces new relevant scal- 
ing fields fl60ll for the fixed points I, EH, and IV, while the fixed 
point II has the same instabilities as in the absence of DDI. 
The new dipolar fixed points in Eqs. ( (38| ) and ( |39"1 > have both 
relevant and irrelevant scaling fields, and the RG flows in the 
vicinity of these fixed points have properties similar to those 
of dipolar ferromagnets with spatial disorder [61 1. We do not 
yet dwell on this point since our analysis so far has neglected 
the relevant single-particle term ( |27] i generated by DDI. Since 
the dipole-dipole interaction does not renormalize to zero in 
the thermal regime, we have to include the single-particle term 
in Eq. ( p7| into our analysis with an a priori unknown cou- 
pling constant h in order to properly investigate the dipolar 
fixed points. 

We analyze the properties of RG flows in the thermal 
regime more carefully in the next section where we consider 
an extensive model for dipolar Bose gases. We also point out 
that the above conclusions hold even if the flow of the chem- 
ical potential is taken into account, i.e., the fixed points and 
their properties remain qualitatively the same. 



VI. COMPLETE RG ANALYSIS AT FINITE 
TEMPERATURES 

We analyze here the properties of dipolar spinor Bose gases 
using the full effective action 5f u ii = S + S + S int , where 
S , S , and 5 int are given by Eqs. ([TJ, ( |2"7] i, and Fur- 
thermore, we allow renormalization of the kinetic-energy term 
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by redefining e k — Z H 2 k 2 /2m with the initial condition 
Z v (0) = 1, In the RG transformation, the kinetic energy 
scales as e k — > e k e lnZ i^~ 2e , and in order to keep the to- 
tal kinetic energy unchanged, we have to rescale fields as 

V*a i/ , Q e?£ ~^ lnZ " W - The quantity lnZ (f) gives rise to 
anomalous dimension 77 which we will discuss in more detail 
once we have the final RG equations at hand. The anomalous 
dimension changes the scaling relations (|9a]i-(|9c|) and gives 
h a nontrivial scaling 



h ^h e z nW, 



(40a) 
(40b) 
(40c) 
(40d) 



where g = g d , g s , c dd . The appearance of nonzero anomalous 
dimension does not change Eq. (|8). 

Since dipole-dipole interactions were found to be rele- 



vant only in the thermal regime (Sect. VB i, we require again 
that the temperature does not flow. We observe that T does 
not acquire any renormalization beyond the trivial rescaling 
even in the case of the augmented action Sf u \\. Assuming 
that In Z (£) does not become too large during the RG flow, 



Eq. ( |40a[ ) renders T to flow to zero. Therefore, both bubble 
diagrams in Fig.[5]give an equal contribution. 

In the presence of the new single-particle operator S , the 
non-interacting Green's function becomes non-diagonal and 
takes the form 



Go, a p(k,u> n ) 



k a kp 



-ihu> n + e k — [i — h k 2 k 2 



(41) 



We note that in the limit h — > 0, Eq. ( pT} reduces to Eq. |6|, 
and if h ^ 0, Go, a p(k, co n ) is non-singular for k — > 0. A 
free propagator analogous to that of Eq. ( pT) has been pre- 
viously considered in the context of dipolar magnets, and 
Eq. ( |4"Tj ) corresponds to the long wave-length limit of the dipo- 
lar propagator of Ref. |52"1 . The difference between the ear- 
lier studies [62] and the present work is the itinerant nature of 
magnetism in spinor Bose gases which gives rise to S only 
through the RG transformation. In the theory of classical mag- 
nets, single-particle interactions similar to S Q represent the ac- 
tual dipole-dipole interaction and phenomenological quartic 
terms are taken to be local interactions 11621 . In particular, we 
see later on that the behavior of h under the RG flow is differ- 
ent between classical dipolar ferromagnets and dipolar spin-1 
Bose gases, even though h has formally the same role in both 
systems. 

Since Go.afi has become nondiagonal, there are new dia- 
grams contributing to the renormalization of coupling con- 
stants. The new diagrams are illustrated in Fig. [9] In order 
to evaluate the angular integrals arising from DDI and nondi- 
agonal Green's function, we compute the RG equations only 
up to the linear order in h . This is a natural approximation, 
since we assume that initially h is small if not vanishing. We 





Figure 9. Additional one-loop diagrams contributing to the RG 
equations due to the nondiagonal part of the Green's function. 



analyze the accuracy of this approximation when we consider 
the fixed points corresponding to the full RG equations. The 
diagrams in Figs.|4j[8} and|9]yield the following contributions: 



dZ„ = 



4m 



c dd (l + 5/i )x ,iAd£, 



IbTth 2 
2 

^ = Y^ C dd( 3 + 5/l o)X ,l Ad ^ 



(42) 
(43) 



= [ - ^(9s + 2s d )(3 + ho) + ^c dd h a ] x ,iA 3 d£, 

(44) 



dg d 



dc dd = 



^ [(3ff s 2 + 6g s g d - S lC 2 d )(l + 2h Q /3) 

+ ( s 29 s c dd + S 3 g d c dd + ,5 4 c^ d )/i ]xo,2 A3 d^, 

(45) 

~[(5g 2 d + 4g d g s + 4.g 2 + 6 5 c 2 dd )(l + 2h /3) 

~ ( S e9d c dd + Vdd)fy)]Xo,2 A3 d^ 

(46) 

1 

_ 2V 2 

-^^d^o]Xo, 2 A 3 d£, (47) 



r2 [(<Vdd + 2 ff d c dd - 4 5s c dd )(l + 2h /3) 



where h = h a A 2 /(e A -n) andx ,„ = l//3(s A -fJ,) n forn = 
1,2. The numerical constants are given by 6 1 = 3047r 2 /45, 
6 2 = 16tt/9, S 3 = 32tt/9, S 4 = 128tt 2 /135, S 5 = 752tt 2 /45, 
S 6 = 2S 3 , S 7 = 3<5 4 , S 8 = 8tt/3, and <5 9 = 5 2 . In the Zeeman 
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basis, the full RG equations take the form 



din Z 

At 
dh o = 

d/x 
df 



i(l 



87TTO 



5h )K 3 A XoA , 
5h )K 3 Ax oa , 



(2-r 1 )»-[-(c 2 +2c )(3 + h ) 



dco 
d^ 



(1 - 2 V )c - [(5cg + 2^ + 



ao c dd^o]-^3 A Xo,i, 
2~ 



(48) 
(49) 

(50) 



a n c; 



dd)(l+ 3^0) 



(a 2 c, 



L dd 



+ "3 c 2 c dd + a4Cdd)^o]" f£: 3 A3 Xo,2. 



(51) 



df 



= (1 - 2 ?/ )c 2 - [(3t| + 6c c 2 + a 5 c§ d )(l + j^ ) 



(a 6 c 2 c dd + a 7 c c dd + a 8 c^ d )/i ]^ 3 A 3 Xo. 



2- 



(52) 



dc. 



del 



df 



(1 



2? ?) c dd " 
X(l 



[(a. 



9 L dd 



2 c c dd 



6c 2 c dd ) 



«io c dd /l o]- pf 3 A3 Xo,2> (53) 



where we have defined the anomalous dimension as 77 = 
d\n.Z v /d£. Comparison with Eqs. (|20), |2T}, ((25), and |26) 
shows that 77 can be thought of as an effective correction to 
the spatial dimension of the system. Numerical constants a { 



are defined as a Q = 87r/3, 



448tt745, a 2 = 32tt/9, 



a 2 /2, 



= 256tt7135, a 5 = 3047r 2 /45, a 6 = a 3 
a 4 /2, a 9 = 



and a 10 = 



a 3 . The RG 



equations (|48)-(|53) are computed up to the linear order ho 
The fixed points corresponding to RG equations d48}-(|53 



are shown in Table III where the dimensionless quantities 
are defined as Hq = A 2 h^/e x , /t* = fi* /e A and c* = 
pe\c*/K 3 A 3 . Table 



III 



shows that fixed points VII and VIII 
correspond to relatively large values of h , and the original 
RG equations (|48)-(f53) are no longer reliable in this region 
since they were calculated only up to the linear order in h . 
Linearized RG equations in the vicinity of fixed points V and 
VI give rise to complex eigenvalues. Similar behavior has 
been previously found in the context of dipolar ferromagnets 
with spatially uncorrected disorder iRTTll as well as in systems 
with long-range-correlated disorder [63 1. We note that the ap- 
pearance of complex eigenvalues for fixed points V and VI 
could be an artifact of our approximations, and therefore we 
concentrate on the properties of RG flows in the vicinity of 
fixed points I-IV where our calculation should capture the es- 
sential physics. 

The Gaussian fixed point I is trivial since it is unstable to 
every direction in the space of the original coupling constants. 
Fixed points II-IV have certain common features such as the 
existence of one marginal scaling field arising from the com- 
bination of h and c dd . This scaling field reflects the existence 
of a continuous line of fixed points for c dd = 0, corresponding 
to an arbitrary value of h$. Since h was originally generated 
by DDI, we have taken Hq — for fixed points with c 



0. 

Fixed points II-IV also have one scaling field directly propor- 



dd 





I II 


111 


IV 


V 


VI 


VII 


VIII 


ho 











-0.355 


-0.357 


1.229 


1.247 


A* 


0.250 


0.250 


0.286 


0.282 


0.279 


0.204 


0.124 


CO 


0.050 


0.109 


0.102 


0.066 


0.070 


0.028 


0.034 




0.088 


-0.031 





0.095 


0.120 


0.046 


0.039 


C dd 











0.007 


-0.009 


-0.002 


0.009 


V* 











-0.012 


0.015 


-0.015 


0.068 



Table III. Dimensionless values of the fixed points corresponding to 
RG equations j48)-{53). All fixed points are unstable and linearized 
RG equations possess complex eigenvalues in the vicinity of fixed 
points V,VI, and VIII. 



tional to c dd . This scaling field is irrelevant for fixed point II 
and relevant for fixed points III and IV. 

To further analyze the behavior of RG flows in the case of 
weak DDI, we simplify the full RG equations by taking both /j 
and h to be critical ll57l . This gives the reduced RG equations 



^ = (1 - 2ry) Co - (5c% + 2cj + a^)^-, 



(54) 



'A 

dc K A 3 

= (1 - 2r;)c 2 - (3 C 2 + 6c c 2 + a 5 c 2 dd ) , (55 ) 

,/ ' ,W _ (1 - 2 »7) c dd - ( a 9 c ld + 2 c c dd + 6c 2 c dd ) K A 



df 



1 ' 



(56) 



where the anomalous dimension is given by 77 = 
j§jjic dd K 3 / f3Ae A . Apart from the anomalous dimension, the 
reduced RG equations are the same as Eqs. d35)-((37). 



The reduced RG equations (|54)-((56) can be used to justify 
the approximation of Sec.|V] where the anomalous dimension 
was neglected altogether. Equations (|54)-(|56) have four fixed 
points I-IV shown in Table|ll|corresponding to vanishing DDI 
(in 3D, we take e = 1). In addition, there are two other non- 
trivial fixed points 



c* a = 0.085, c 2 = 0.114, c* dd = 0.015, (57) 
= 0.150, c dd = -0.008. (58) 



c* a = 0.090 



Comparing Eqs. ((38) and ((39) to Eqs. ((57) and ((58), one ob- 
serves that the effect of anomalous dimension is negligible. 
Furthermore, Table [ill] shows that the values of r\ at fixed 
points are small compared to unity and hence the anomalous 
dimension has only a small effect on the fixed point structure 
of the full RG equations (|48)-(|53). 

In the vicinity of fixed points I-IV, linearized RG equations 
again give rise to one scaling field that is directly proportional 
to the dipole-dipole interaction. In the case of fixed points I, 
III, and IV, this scaling field is relevant and DDI introduces 
an additional instability. To quantify this instability, we de- 
fine crossover exponents [46 1 <f> i = Aj/A DDI , where A DDI is 
the eigenvalue corresponding to the DDI-induced scaling field 
and Aj 2 are the two remaining eigenvalues corresponding to 
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the flows depicted in Fig. [7] The crossover exponent <p i in- 
dicates the relative importance of a given scaling field s i with 
respect to the DDI induced instability ||60ll . When the absolute 
value of the crossover exponent becomes smaller than unity, 
the DDI-dependent instability dominates the RG flow with re- 
spect to Sj. 





I 


11 


III 


IV 


01 


1 


5.498 


-1.069 


-1.667 


4>2 


1 


-8.464 


0.235 


-0.333 



Table IV. Crossover exponents corresponding to fixed points I-IV in 
Table [TT| For the fixed point II, we have replaced A DDI with |A DDI 
since the DDI scaling field is irrelevant at fixed point II. 



The crossover exponents are shown in Table IV and they 
indicate that the properties of fixed point II are largely unaf- 
fected by the dipole-dipole interaction. On the other hand, 
fixed points I, II, and IV are susceptible to the runaway flows 
induced by DDI. To explore these runaway flows, we integrate 
the RG equations numerically in the vicinity of fixed points I- 
IV. We take Cdd(0) to be small and positive, which is the phys- 
ically relevant case. We find that the system exhibits the two 
runaway flows discussed in Sect. IV (see Fig[7). However, the 
mean-field regime governed by the fixed point IV becomes 
unstable and RG flows starting from this region correspond to 
the ferromagnetic runaway flow. The dipole-dipole interac- 
tion tends to increase under the ferromagnetic runaway flow, 
whereas under the antiferromagnetic runaway flow DDI renor- 



malizes to zero. This is demonstrated in Figs. 10 and 1 1 where 
we show the sign of c 2 (Fig. 10 1 and the magnitude of c, 



da 



(Fig. 1 1 1 in the asymptotic limit of RG flows. For each initial 



point (c , c 2 , c dd ), we integrate the RG equations up to value 
£ c , where £ c is given by the condition |c 2 (^ c )| = 10- Since 
both c and c 2 diverge along the runaway flows, the precise 
value of |c 2 (^ c )| which determines l c is unimportant. It only 
needs to be large enough to illustrate the general tendency as- 
sociated with the two runaway flows: for AFM runaway flow 
DDI renormalizes to zero whereas for FM runaway flow DDI 
tends to grow. 



Figures 10 and 11 demonstrate the existence of a phase 
transition separating an antiferromagnetic spin singlet con- 
densate and a ferromagnetic condensate with anisotropic 
dipole-dipole interactions and suppressed fluctuations in the 
magnitude of local spin (dipolar ferromagnetic condensate). 
The dipole-dipole interaction favors spatial modulation in the 



local magnetization and in contrast to Sect. IV B the equilib- 
rium phase corresponding to the ferromagnetic runaway does 
not feature all atoms in the same hyperfine spin state. Fur- 
thermore, even weak DDI renders antiferromagnetic and fer- 
romagnetic condensates unstable towards formation of dipolar 
ferromagnetic condensate. 

To verify that the anisotropic kinetic energy term repre- 
sented by h does not change the previous conclusions, we lift 
the assumption h Q = and integrate RG equations (|48|l-(|53) 
in the critical region corresponding /i = 0. We take h ini- 
tially zero since the original action does not contain the term 
in Eq. ( |27] ). We obtain RG flows identical to those depicted in 




Figure 10. (Color online) Sign of the dimensionless coupling con- 
stant c 2 (£ c ) for each initial point (c ,c 2 ,c dd ). The dimensionless 
coupling constants are defined as = I3e\cil K D A D ', i — 0, 2, dd. 
Reduced RG equations (|54}-l|56[) are integrated up to l c , where £ c is 
determined by the condition |c 2 (£ c )| = 10. The initial values of c dd 
are (a) c dd (0) = 0.01, (b) c dd (0) = 0.025, (c) c dd (0) = 0.05, and 
(d) c dd (0) = 0.1. 
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Figure 1 1 . (Color online) Value of the dimensionless coupling con- 
stant c dd (£ c ) for different initial points (c ,c 2 ,c dd ). The reduced 
RG equations are integrated up to £ c , where £ c is determined by the 



condition c 2 ( 



= 10. The initial values of c dd are the same as 



in Fig. [10] namely (a) c dd (0) 
c dd (0) = 0.05, and (d) c dd (0) = 



= 0.01, (b) c dd (0) 
0.1. 



0.025, (c) 



Fig.(10jand 
we find that 
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For the physically relevant case c dd (0) > 0, 
7Iq grows under both runaway flows. Under the 
antiferromagnetic runaway flow, the growth of h is somewhat 
slower than in the case of ferromagnetic runaway flow. 

We conclude this section by analysing the stability of the 
two additional fixed points given in Eqs. ( |57| ) and ((58). As 
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discussed in Sec. VB the RG flow in the vicinity of these 
new fixed points has somewhat unconventional properties that 
are similar to those discussed in Refs ||6T1 |63ll . When the 
RG equations (|54"|)-(|56l) are linearized in the vicinity of fixed 
points, the eigenvalues of the resulting matrix consist of one 
real eigenvalue and a pair of complex-conjugated eigenvalues 
in the case of physically relevant fixed point in Eq. ( f57] i. The 
real parts of eigenvalues determine the stability of the fixed 
point [63], and we find that the fixed point in Eq. ( f57"| is un- 
stable and practically unattainable since the complex eigen- 
values have positive real parts. Furthermore, RG flows near 
the fixed point in Eq. ( |57| ) are not markedly different from the 
RG flows near fixed points I-IV. We find that RG flows start- 
ing from positive c dd remain positive and since the bare value 
of DDI coupling c dd (0) = i^Q^g 2 / An is positive, the fixed 
point in Eq. ( |58] l corresponding to negative DDI coupling is 
not relevant for physical systems. For completeness, we note 
that linearized RG equations for the fixed point in Eq. ( |58) l 
give rise to real eigenvalues, two of which are positive. 



VII. EFFECTS ASSOCIATED WITH THE ANISOTROPIC 
DISPERSION 

To analyze the properties of the anisotropic kinetic energy 
term in Eq. ( |27| >, we restrict to the noninteracting limit and 
consider only the single-particle part So+S' Q . Since the effects 
associated with dipole-dipole interactions become important 
only at finite temperatures, we neglect all nonzero Matsubara 
frequencies in the correlation function Go, Q /3 = ( , 4>a' t l>p) an d 
obtain 



1 



/' 



k 2 ' 



1 



/U — h k 2 k 2 



(59) 

The anisotropic dispersion does not affect local spin order, 
but can nevertheless give rise to nematic order described by a 
tensor order parameter 



1 



(60) 



Since we consider a uniform system with constant total den- 
sity, our definition of Q s is analogous to that of Ref. Il64l . 

In the non-interacting limit we have (Q s (k)) = G a p(k). 
The nematic order is associated with the eigenvalues and 
eigenvectors of (Q s ). Following Ref. 1641 . we define the ne- 
matic director n to be the eigenvector corresponding to the 
largest eigenvalue of (Q s )- Since (Q s ) has to be positive 
semidefinite, we impose a condition \h \ < h 2 /2m. This 
condition is also physically motivated since initially in the 
RG calculation, we assumed h to be small compared to ki- 
netic energy. Without losing generality we may also take 
/j = 0. For h > we obtain n(k) = k, corresponding 
to a hedgehog (monopole) in the momentum space. In the 
case h < 0, there are two linearly independent nematic di- 
rectors which can be taken to be n 1 (k) = (—k z ,0, k x )/k xz 
and n 2 (fe) = (-k y ,k x ,0)/k xy , where k xa = ^Jkl + k 2 a . 
Such nematic directors correspond to vortices in momentum 
space. 



One of the experimental manifestations of dipolar interac- 
tions in spin-1 Bose gases is periodic spatial modulation in the 
local magnetization [27]. Since we have considered only the 
noninteracting case in this section, the local spin order van- 
ishes and our analysis cannot be directly compared with the 
experiment of Ref. 11271 . To study the nematic order associ- 
ated with h in real space, we Fourier transform Go, a/3- The 
nematic director n(r) corresponds to the hedgehog for h < 
and the two vortex solutions for h > 0. Note that the hedge- 
hog in real space corresponds to a vortex in momentum space 
and vice versa. The nematic ordering and the corresponding 
textures can in principle be measured in real space utilizing 
the atomic birefringence [65 66). 



VIII. DISCUSSION 

In this work, we have analyzed the properties of dipo- 
lar spin-1 Bose gases using momentum shell renormalization 
group and taking into account all one-loop diagrams. In the 
absence of magnetic dipole-dipole interactions, our RG anal- 
ysis complements the previous RG studies of low-dimensional 
spinor Bose gases at T = ll39l as well as the phenomeno- 
logical analysis of Ref. 11381 . Similarly to Ref. (39], we found 
two runaway flows corresponding to the formation of either a 
condensate of spin singlet pairs or a fully spin-polarized scalar 
condensate. The absence of stable fixed points in the regions 
corresponding to runaway flows is sometimes a manifestation 
of a first-order transition l54l 1531 . In the case of antiferro- 
magnetic runaway flow, it would be interesting to study if sta- 
ble fixed points arise when an additional field representing the 
spin singlet pairs is introduced [39]. We believe that the pos- 
sible first-order transition associated with the ferromagnetic 
runaway flow could be akin to the fluctuation-induced first- 
order transition in type I superconductors ll67l . For spinor 
Bose gases, the analogue of an intrinsic fluctuating magnetic 
field is given by the fluctuating Berry phase associated with 
the local magnetization J68). 

In the zero-temperature limit, we found that the dipole- 
dipole interaction renormalizes to zero and does not induce 
any additional instabilities. At finite temperatures, we ana- 
lyzed the limit where thermal fluctuations dominate quantum 
fluctuations. We found that the pair condensate is unaffected 
by the dipole-dipole interactions which eventually renormal- 
ize to zero. On the other hand, both antiferromagnetic and fer- 
romagnetic condensates become unstable and the system ex- 
hibits an instability similar to the ferromagnetic runaway flow 
in the absence of DDI. In principle the magnetic dipole-dipole 
interaction can be transformed to an external vector potential 
such that the local spin of the gas S couples linearly to the curl 
of the vector potential [48]. This transformation gives rise to 
an alternative RG scheme which could provide further insight 
to the runaway flow induced by the DDI and its connection 
to the potential first-order transition. Also, the role of higher 
order terms beyond the one-loop approximation should be ex- 
plored and one possible route to accomplish this task could be 
the functional renormalization group, from which the current 
RG equations arise in principle as an approximation Il34ll36l . 
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Since the lifetime of ultracold atomic gases is limited and 
the spin-spin interactions are relatively weak, it is not clear 
to what extent the current experiments |26l l27ll are able to 
explore the true thermal equilibrium of the system. How- 
ever, even if the experimentally attainable physics of spinor 
Bose gases eventually turns out to be inherently out of equilib- 
rium, understanding the corresponding equilibrium systems is 
still a prerequisite for the exploration of the non-equilibrium 
situation. The experimentally relevant atomic species 23 Na 
and 87 Rb give rise to bare coupling constants that belong to 
the regions of antiferromagnetic and ferromagnetic conden- 
sates, respectively. In the absence of dipole-dipole interac- 
tions, the critical properties of both of these condensates are 
determined by the S'[/(3)-symmetric fixed point discussed in 
Section IV B When dipole-dipole interactions are taken into 
account, ferromagnetic and antiferromagnetic condensates be- 
come in principle unstable and the true equilibrium is deter- 
mined by the ferromagnetic runaway flow. However, the life- 
time of atomic gases can limit the possibilities of observing 



this crossover from the critical behavior determined by the 
non-dipolar fixed points of Sec. |IV B| and [VI] to the thermo- 
dynamic equilibrium determined by DDI. 

Recently observed optical Feshbach resonances iPffl as well 
as the proposed microwave induced Feshbach resonances ll42l 
provide in principle means to fully explore the phase diagrams 
studied in this work. Alternatively, the phase diagrams in the 
absence of DDI could also be realized in the molecular super- 
fluid phase of p-wave resonant Bose gases ll69l . 



ACKNOWLEDGMENTS 

We would like to thank E. Demler, D. Podolsky, and 
P. Kuopanportti for insightful comments. This work was sup- 
ported by the Academy of Finland (VP and MM), Emil Aalto- 
nen Foundation (VP and MM), and Harvard-MIT CUA (VP). 



[1] A. Micheli, G. K. Brennen, and P. Zoller, Nat. Phys. 2, 341 
(2006). 

[2] T. Donner, S. Ritter, T. Bourdel, A. Ottl, M. Kohl, and 

T. Esslinger, Science 315, 1556 (2007). 
[3] R. Jordens, N. Strohmaier, K. Giinter, H. Moritz, and 

T. Esslinger, Nature 455, 204 (2008). 
[4] G.-B. Jo, Y.-R. Lee, J.-H. Choi, C. A. Christensen, T. H. Kim, 

J. H. Thywissen, D. E. Pritchard, and W. Ketterle, Science 325, 

5947 (2009). 

[5] A. Gorshkov, M. Hermele, V. Gurarie, C. Xu, P. S. Julienne, 

J. Ye, P. Zoller, E. Demler, M. D. Lukin, and A. M. Rey, 

Nat. Phys. 6, 289 (2010). 
[6] T.-L. Ho, Phys. Rev. Lett. 81, 742 (1998). 
[7] T. Ohmi and K. Machida, J. Phys. Soc. Jpn. 67, 1822 (1998). 
[8] U. Leonhardt and G. E. Volovik, JETP. Lett. 72, 66 (2000). 
[9] H. T. C. Stoof, E. Vliegen, and U. Al Khawaja, Phys. Rev. Lett. 

87, 120407 (2001). 
[10] U. Al Khawaja and H. Stoof, Nature 411, 918 (2001). 
[11] C. M. Savage and J. Ruostekoski, Phys. Rev. A 68, 043604 

(2003). 

[12] M. Takahashi, S. Ghosh, T. Mizushima, and K. Machida, 

Phys. Rev. Lett. 98, 260403 (2007). 
[13] Y. Kawaguchi, M. Nitta, and M. Ueda, Phys. Rev. Lett. 100, 

180403 (2008). 

[14] V. Pietila' and M. Mottonen, Phys. Rev. Lett. 103, 030401 
(2009). 

[15] J. A. M. Huhtamaki, M. Takahashi, T. P. Simula, T. Mizushima, 
and K. Machida, Phys. Rev. A 81, 063623 (2010). 

[16] J. A. M. Huhtamaki and P. Kuopanportti, Phys. Rev. A 82, 
053616(2010). 

[17] S. Mukerjee, C. Xu, and J. E. Moore, Phys. Rev. Lett. 97, 
120406 (2006). 

[18] V. Pietila, T. P. Simula, and M. Mottonen, Phys. Rev. A 81, 
033616 (2010). 

[19] J. Stenger, S. Inouye, D. M. Stamper-Kurn, H.-J. Miesner, A. P. 

Chikkatur, and W. Ketterle, Nature 396, 345 (1998). 
[20] H. Schmaljohann, M. Erhard, J. Kronjager, M. Kottke, S. van 

Staa, L. Cacciapuoti, J. J. Arlt, K. Bongs, and K. Sengstock, 

Phys. Rev. Lett. 92, 040402 (2004). 



[21] M.-S. Chang, C. D. Hamley, M. D. Barrett, J. A. Sauer, K. M. 
Fortier, W. Zhang, L. You, and M. S. Chapman, Phys. Rev. Lett. 
92, 140403 (2004). 

[22] L. E. Sadler, J. M. Higbie, S. R. Leslie, M. Vengalattore, and 
D. M. Stamper-Kurn, Nature 443, 312 (2006). 

[23] A. Griesmaier, J. Werner, S. Hensler, J. Stuhler, and T. Pfau, 
Phys. Rev. Lett. 94, 160401 (2005). 

[24] T. Lahaye, T. Koch, B. Frolich, M. Fattori, J. Metz, A. Gries- 
maier, S. Giovanazzi, and T. Pfau, Nature 448, 672 (2007). 

[25] T. Koch, T. Lahaye, J. Metz, B. Frohlich, A. Griesmaier, and 
T. Pfau, Nat. Phys. 4, 218 (2008). 

[26] M. Vengalattore, S. R. Leslie, J. Guzman, and D. M. Stamper- 
Kurn, Phys. Rev. Lett. 100, 170403 (2008). 

[27] M. Vengalattore, J. Guzman, S. R. Leslie, F Serwane, and 
D. M. Stamper-Kurn, Phys. Rev. A 81, 053612 (2010). 

[28] K. G. Wilson and J. Kogut, Phys. Rep. 12, 75 (1974). 

[29] R. Shankar, Rev. Mod. Phys. 66, 129 (1994). 

[30] D. S. Fisher and P. C. Hohenberg, Phys. Rev. B 37, 4936 (1988). 

[31] E. B. Kolomeisky and J. P. Straley, Phys. Rev. B 46, 11749 
(1992). 

[32] E. B. Kolomeisky and J. P. Straley, Phys. Rev. B 46, 13942 
(1992). 

[33] M. Bijlsma and H. T. C. Stoof, Phys. Rev. A 54, 5085 (1996). 
[34] C. Wetterich, Nucl. Phys. B 352, 529 (1991). 
[35] J. O. Andersen and M. Strickland, Phys. Rev. A 60, 1442 
(1999). 

[36] J. O. Andersen, Rev. Mod. Phys 76, 599 (2004). 

[37] S. Diehl, S. Floerchinger, H. Gies, J. M. Pawlowski, and 
C. Wetterich, Ann. Phys. (Berlin) 522, 615 (2010). 

[38] K. Yang, "Phase Diagrams of Spinor Bose Gases," (unpub- 
lished), arXiv:0907.4739, (2009). 

[39] A. K. Kolezhuk, Phys. Rev. A 81, 013601 (2010). 

[40] S. S. Natu and E. J. Mueller, "Pairing, Ferromagnetism, and 
Condensation of a normal spin-1 Bose gas," (unpublished), 
arXiv:1101.5639, (2011). 

[41] C. D. Hamley, E. M. Bookjans, G. Behin-Aein, P. Ahmadi, and 
M. S. Chapman, Phys. Rev. A 79, 023401 (2009). 

[42] D. J. Papoular, G. V. Shlyapnikov, and J. Dalibard, Phys. Rev. 
A81,041603(R)(2010). 



13 



[43] R. W. Cherng and E. Dernier, Phys. Rev. Lett. 103, 185301 

(2009) . 

[44] T. Lahaye, C. Menotti, L. Santos, M. Lewenstein, and T. Pfau, 

Rep. Prog. Phys. 72, 126401 (2009). 
[45] G. E. Astrakharchik and Y. E. Lozovik, Phys. Rev. A 77, 

013404 (2008). 

[46] S.-K. Ma, Modern Theory of Critical Phenomena (Addison- 

Wesley, Reading, MA, 1976). 
[47] D. W. Snoke, Solid State Physics: Essential Concepts 

(Addison-Wesley, Boston, 2008). 
[48] D. Belitz and T. R. Kirkpatrick, Phys. Rev. B 81, 184419 

(2010) . 

[49] S. Yi and H. Pu, Phys. Rev. Lett. 97, 020401 (2006). 
[50] Y. Kawaguchi, H. Saito, and M. Ueda, Phys. Rev. Lett. 97, 
130404 (2006). 

[51] J. P. Burke, C. H. Greene, and J. L. Bohn, Phys. Rev. Lett. 81, 
3355 (1998). 

[52] H. Bruus and K. Flensberg, Many-Body Quantum Theory in 
Condensed Matter Physics (Oxford University Press, Oxford, 
2004). 

[53] S. Sachdev, Quantum Phase Transitions (Cambridge University 
Press, Cambridge, 1999). 

[54] J. Rudnick, Phys. Rev. B 18, 1406 (1978). 

[55] D. J. Amit and V. Martin-Mayor, Field Theory, the Renormal- 
ization Group, and Critical Phenomena, 3 Td Ed. (World Scien- 
tific, Singapore, 2005). 

[56] T.-L. Ho and S. K. Yip, Phys. Rev. Lett. 84, 4031 (2000). 

[57] We take mass terms such as n to be critical (zero in the case of 
jj,) since they correspond to a temperature instability that tends 
to drive the system away from the critical point under the RG 
transformations, see Ref. 15 81 1631 . In practice, this means that 
temperature has to be tuned in order to reach the critical point. 



When the coupling constant ho arising from an anisotropic dis- 
persion [see Eq. l |27[ (] is present, we have two parameters asso- 
ciated with flows that drive the system away from the critical 
region under successive applications of the RG transformation. 
Hence the critical points examined in Section [VT] are examples 
of bicritical points which can be reached by tuning both p. and 
h . Biciritical fixed points occur also, e.g., in anisotropic mag- 
netic materials F70l . 

[58] A. Aharony in Phase Transitions and Critical Phenomena, 
edited by C. Domb and M. S. Green (Academic Press, New 
York, 1976), Vol. 6. 

[59] E. Eisenberg and E. H. Lieb, Phys.Rev. Lett. 89, 220403 (2002). 

[60] The scaling fields are defined as projections of the original cou- 
pling constants (c , c^, c dd ) on the left-eigenvectors of the ma- 
trix corresponding to the linearized RG equations. 

[61] A. Aharony, Phys. Rev. B 12, 1049 (1975). 

[62] A. Aharony and M. E. Fisher, Phys. Rev. B 8, 3323 (1973), and 
references therein. 

[63] A. Weinrib and B. I. Halperin, Phys. Rev. B 27, 413 (1983). 

[64] E. J. Mueller, Phys. Rev. A 69, 033606 (2004). 

[65] I. Carusotto and E. J. Mueller, J. Phys. B.: At. Mol. Opt. Phys. 
37, SI 15 (2004). 

[66] J. M. Higbie, L. E. Sadler, S. Inouye, A. P. Chikkatur, S. R. 

Leslie, K. L. Moore, V. Savalli, and D. M. Stamper-Kurn, 

Phys. Rev. Lett. 95, 050401 (2005). 
[67] B. I. Halperin, T. C. Lubensky, and S.-K. Ma, Phys. Rev. Lett. 

32, 292 (1974). 

[68] T.-L. Ho and V. B. Shenoy, Phys. Rev. Lett. 77, 2595 (1996). 
[69] L. Radzihovsky and S. Choi, Phys. Rev. Lett. 103, 095302 
(2009). 

[70] P. M. Chaikin and T. C. Lubensky, Principles of condensed mat- 
ter physics (Cambridge University Press, Cambridge, 2000). 



